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Abstract 

A non-linear theory for the plastic deformation of prismatic bodies is 
constructed which interpolates between Prandtl's linear soap-film approx- 
imation and Nadai's sand-pile model . Geometrically Prandtl's soap film 
and Nadai's wavefront are unified into a single smooth surface of constant 
mean curvature in three-dimensional Minkowski spacetime. 
Keywordsrtorsion, prismatic body, sandpile analogy, Minkowski space- 
time. 
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1 Introduction 

Prandtl's soap film Q] and Nadai's [2] sand heap analogies remain important 
tools for analysing the torsional loads and plastic deformations of cylindrical 
shafts [H in this paper we propose a generalisation of Prandtl's model 
in which the eponymous potential <p(x, y) he introduced no-longer satisfies the 
linear Laplace equation but rather a non-linear equation known as the Born- 
Infeld equation. This has the property that the Prandtl potential 4>{x, y) is 
continuous but nevertheless the torsional stress |V</>| never exceeds the plastic 
limit. The model also permits a remarkable unified continuous geometrical 
analogue combining both the soap-film and sand pile analogies as limiting cases. 
According to this unified analogy 4>(x, y) defines a spacelike maximal surface in 
an auxiliary 2+1 dimensional spacetime. 
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2 The Prandtl-Nadai Construction 



In what follows we provide, for completeness, a brief resume of Prandtl's soap- 
film and Nadia's sand heap analogy In the absence of body forces, the governing 
equations are 

diTij = (1) 

where = Tji is the stress tensor (often often denoted in Engineering texts 
<7jj or Tij cf. [3]). The strain tensor is defined in terms of displacements Ui 
as eij = d{Uj + djUi. In linear theory for an isotropic substance Hooke's Law 
becomes 

«- ^^-1%^. (2) 

where v is Poisson's ratio and E is Young's modulus. Following Saint- Venant 5. 
we assume that the material is isotropic and confined to a prismatic cylinder 
with base B whose generators are parallel to the z axis and orthogonal to 
the plane curve 7 = dB . The displacements (u,v,w) are assumed to satisfy 
u = —Oyzv = +9xzw — w(x,y), where 9 is a constant. The function w(x,y) 
is called the warp function and it gives the displacement parallel to the axis. 
These assumptions imply 

&xx — &yy — e zz = &xy = O7 &yz = dyW + Ox, e zx = d x W 9y. (3) 

Hooke's law implies that the diagonal components of the stress tensor (the 
elongations) vanish and therefore equation ([1} gives 

d x T xz + d x T yz = (4) 

(this only requires d z T zz — ). We can now introduce the Prandtl potential 

T xz = d y (j) , T yz = -d x (/) . (5) 

We have from (J5J) di(f)Ti Z = ,i = x, y which implies that the component of the 
force in the z direction on the level sets of the Prandtl potential, whose normal 
is n^jj , vanishes. Thus any level set may act as a free boundary. It is convenient 
to set (j> — on the boundary 7 = dB. The total torque is given by 

T = I (xT yz - yT xz ) dxdy = 2 4>dxdy (6) 

J B JB 

on the use of the divergence theorem and the boundary conditions. If G = 
E/2(l + v), and we use Hooke's law we get 

T yz = G{9y + d y w) , T xz = G{-6x + d x w) . (7) 

Thus the Prandtl potential satisfies Poisson's equation 

dx 2 dy' 2 ^ ^ 
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This must be solved in B, subject to the boundary condition that </> = on the 
boundary 7 = dB. Of course Poisson's equation arises in electrostatics where 
the electric field, and electrostatic potential (j> satisfy E = — V , V ■ E = p, 
where p is the density of electric charge. In the case of Prandtl's equation 
© we have p = 2G9 which corresponds to constant charge density. Prandtl 
adopted a different analogy. He regarded <j> — 4>{x,y) as a height function in 
some auxiliary three dimensional space Euclidean space E 3 with coordinates 
(xi, X2, Xs) — (x, y, 4>) then the surface is approximately of constant mean cur- 
vature such as would be adopted by a soap-film with a constant pressure differ- 
ence. The torque is then proportional to the volume between the approximate 
soap film and the base B. Note that if the analogy were exact then the Prandtl 
potential would satisfy the Young-Laplace equation 

dx M +d M = _ 2Ge (9) 
VT+WW 2 Vi + W 

2.1 Non-linear theory 

If the total shear stress exceeds the elastic limit k, then the material deforms 
plasticly. Thus a plastic upper bound must be satisfied 

t 2 z + t 2 z = \\7^\ 2 <k\ (10) 

where k is the plastic limit and in the plastic region the Prandtl potential is 
assumed to satisfy the Eikonal equation 

|V0| 2 = fc 2 . (11) 

According to Nddai 's sandpile analogy [2J [5J [7] one may regard a Prandtl poten- 
tial satisfying the Eikonal equation (TTT|) as the height function <fi of a sand-heap 
or sand-pile of constant angle of repose a = tan^ 1 k located in same Euclidean 
space E 3 . The strategy for solving the for 4> adopted by Prandtl is to erect over 
the base B, a sand-pile, that is the solution </>Eikonai of the Eikonal equation 
(TrT|) . One then also erects the approximate soap film, that is the graph of the 
solution of Poisson's equation 0p O isson- The solution adopted by Prandtl 
is then to take the minimum of 4>p isson and 4>Eikonai- In other words, linear 
theory is assumed valid under the tent. 

3 Spacetime Interpretation 

We begin by noting that the the Eikonal equation (fTTj) may be regarded as the 
defining a wave-front, null or characteristic surf ace N : u(x,y,t) — where 

u = t-4>(x,y)=0 (12) 

in 3- dimensional Minkowski spacetime E 2,1 with coordinates x,y,t and metric 

ds 2 =dx 2 +dy 2 -c 2 dt 2 , (13) 
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where c = 1/fc is the velocity of light. In other words, 



- ^{d t uf + {d x uf + [d v uf = . (14) 

In general, the plastic upper bound (jlOl) must hold. However the transition be- 
tween linear and plastic behaviour is observed to be smoother than the abrupt 
change envisaged in Prandtl's theory. Now (fit))) implies that the surface E C E 2 ' 1 
has a normal which is non-spacelike. If inequality holds in (|10p then S given by 
(|12p is a spacelike surface. Thus one considers a smooth surface whose normal 
may become null, i.e., may satisfy (lll[) in a continuous fashion. It is natural 
therefore to replace the linear Poisson equation with a non-linear equation which 
interpolates between the Poisson equation and the Eikonal equation and which 
has an interpretation more in keeping with the Minkowski spacetime frame- 
work. On might try replacing Poisson's equation by the Young-Laplace equa- 
tion ([9]) but this does not fit well with our interpretation the Eikonal equation. 
A geometrically better motivated suggestion is to postulate the relativistically 
covariant analogue of the Young-Laplace equation: 

d x - dx<t> +d y - = -2G6, (15) 

which arises as the Euler- Lagrange equations of the functional 



y/i - c 2 |V</>| 2 + 2Gc 2 <9</> dxdy . (16) 



The first term in (ITBl is the area of any spacelike surface E with edge 7 = dB 
and the second the spacetime volume between E and the spacelike hyperplane 
of constant time t = 0. It follows that (|15|) describes a spacelike surface of 
constant mean curvature in three dimensional Minkowski spacetime E 2 ' 1 . The 
difference in sign in the arguments inside the square roots in the denominators 
of ([9]) and (ITSj) is intended to enforce the plastic bound (|TU| . 



4 Born and Infeld's non-linear Electrodynamics 

Just as Poisson's equation (JSJ has an electrostatic analogue, so does our pro- 
posed replacement (fT5l) . The equation (fT5t and associated variational prin- 
ciple with functional (| 16[) admit a similar physical interpretation [5]. Born 
and later Born and Infeld [HJ [TO] suggested a non-linear version of electrody- 
namics: a modification of Maxwell's equations in which there is a maximum 
electric field strength. In the electrostatic situation the electric field strength 
E electric induction D and charge density p satisfy E = — V0V • D = pD = 
E/ — E 2 /6 2 ). In two dimensions, these equation, with p = 2G9, b = 1/c, 
coincide with (fl5|) and were first studied in the context of Born-Infeld electro- 
statics by Pryce [TT1H2]. 
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Figure 1: Comparison between the two models. Solid line: Prandtl-Nadai the- 
ory, dashed line d(f>/dr = GOrk/y/k 2 + G 2 9 2 r 2 : our model. Maximal deviation 
at elasto-plastic boundary. The dashed curve also yields the constitutive law if 
we subsitute 7 = rO for the strains. 



5 Comparison in the circular case 

If the cross section is circular, the standard case has 

GO 2 , „ _ _ k 



-r 



0<r<— (17) 



Y 2 YU , _ _ G0 

= -*r+^ + A>. r>±. (18) 

The inner region is the Prandtl regime. The outer region is the Nadai's sand 
pile regime. The constant c/)q is chosen so that <j>(r = R) = 0, where R is the 
radius of the cylinder. Thus if R < then <fr — ^-R 2 and Prandtl's solution 

is 4> = -\GB{r 2 - R 2 ) while if R > then cf> = kR— and the sandpile 
solution is conical with <j> = —k(r — R) . In the Born-Infeld case one has a single 
unified formula 



, k 2 ( / G 2 9 2 R 2 / G 2 2 r 2 \ 

For small r, we expand both square roots and recover Prandtl's solution at 
lowest order. For large R and hence large r we ignore the one inside the square 
roots and the solution approachs the conical sand pile solution. Note that while 
4>{R) = in both solutions, in general the values of 0(0) are not the same. 
The predictions of both theories are compared in Figure [U The total torque 
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T = 4-7T r0 dr is given by 



fc 2 (\ a , G 2 8 2 R\i k 2 , , G 2 9 2 2 ,,,n 

6 Interpretation in Minkowski spacetime 

We have 

(^>* -«■-»*-;£■ < 21 » 

Thus if we think of (0, x, y) as coordinates for three dimensional (x, y, 0) Minkowski 
spacetime with metric given in (|13|). then should be thought of as the time 
coordinate. Our solutions is now seen as a hyperboloid of constant spacetime 
distance from the point (<j>o, 0, 0), i.e the analogue in Minkowski space of a sphere 
in Euclidean space. To get from the sphere to the hyperbola one may "Wick 
Rotate", i.e., set 

c(cf> - <j> ) = i(z - z )G6/k = io- 1 , (22) 

so that (z — zq) 2 + x 2 + y 2 — a 2 . The fact that the mean curvature (in the 
Lorentzian sense) is constant is now obvious because the surface is invariant 
under Lorentz transformations about (</>o,0, 0) and the mean curvature is a 
Lorentz scalar. 
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